UNCLASSIFIED 


_ AD  NUMBER _ 

AD872813 

LIMITATION  CHANGES 
TO: 

Approved  for  public  release;  distribution  is 
unlimited. 


FROM: 

Distribution  authorized  to  U.S.  Gov't,  agencies 
and  their  contractors;  Critical  Technology;  JUL 
1970.  Other  requests  shall  be  referred  to  U.S. 
Army  Aviation  Materiel  Laboratories  ,  Fort 
Eustis  ,  Virginia  23604.  This  document  contains 
export-controlled  technical  data. 


_ AUTHORITY 

USAAMRDL  ltr,  23  Jun  1971 


THIS  PAGE  IS  UNCLASSIFIED 


u 


Cv> 

IN 

00 

A 


HSAAVLABS  TECHNICAL  REPORT  69-82 

NONLINEAR  FREE  VIBRATIONS  OF  THIN,  CIRCULAR 
CYLINDRICAL  SHELLS 


>- 

c 


o 

o 

1.  Mayers 

LiJ 

1.  6.  Wren 

■ 

cr> 

U— 

- 

s 

a 

c 

C  1 

CZJ 

July  1970 

U.  S.  ARMY  AVIATION  MATERIEL  LABORATORIES 
FORT  EOSTIS,  VIRGINIA 


CONTRACT  DA  44-177-AMC-115(T) 
DEPARTMENT  OF  AERONAUTICS  AND  ASTRONAUTICS 
STANFORD  UNIVERSITY 
STANFORD,  CALIFORNIA 


This  document  it  subject  to  special 
export  controls,  and  each  transmittal 
to  foreign  governments  or  foreign 
nationals  may  be  made  only  with 
prior  approval  of  U.  S.  Army  Aviation 
Materiel  Laboratories,  Fort  Eustit. 
Virginia  23604. 


Disclaimers 


The  findings  in  this  report  are  not  to  be  construed  as  an  official  Depart¬ 
ment  of  the  Army  position  unless  so  designated  by  other  authorized 
documents. 

When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  related  Government 
procurement  operation,  the  United  States  Government  thereby  incurs  no 
responsibility  nor  any  obligation  whatsoever;  and  the  fact  that  the 
Government  may  have  formulated,  furnished,  or  in  any  way  supplied  the 
said  drawings,  specifications,  or  other  data  is  not  to  be  regarded  by 
implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights  or  permission,  to 
manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way 
be  related  thereto. 


Disposition  Instructions 

Destroy  this  report  when  no  longer  needed.  Do  not  return  it  to  the 
originator. 


I 


DEPARTMENT  OF  THE  ARMY 

HCAOQUARTCnS  US  ARMY  AVIATION  MATBML  LABONATOMCS 
PONT  CUSTtS.  VINO  IN  IA  1M04 


This  program  was  carried  out  under  Contract  DA  44- 177-AMC- 115(T)  with 
Stanford  University. 

The  data  contained  in  this  report  are  the  result  of  research  conducted 
to  evaluate  the  influence  of  higher  order  and  nonlinear  effects  on  the 
free  vibration  behavior  of  thin,  circular  cylindrical  sheila. 

The  report  has  been  reviewed  by  the  U.t*.  Army  Aviation  Material  Labora¬ 
tories  and  is  considered  to  be  technically  sound.  It  la  published  for 
the  exchange  of  information  and  the  stimulation  of  future  research. 


Task  IF162204A17002 
Contract  DA  M»-177-AMC-113(T) 
U8AAVLAB8  Technical  Report  69-62 

July  1970 


NONLINEAR  FREE  VIBRATIONS  OF  THIN,  CIRCULAR 
CYLINDRICAL  SHELLS 

by 

J.  Meyers  and  B.  0.  Wrenn 


Prepared  by 
Stanford  University 

Department  of  Aeronautics  and  Astronautics 
Stanford,  California 


for 

U.  S.  ARMY  AVIATION  MATERIEL  LABORATORIES 
FORT  EUSTI8,  VIRGINIA 


This  document  is  subject  to  special  export  controls,  and  each 
transmittal  to  foreign  governments  or  foreign  nationals 
may  be  made  only  with  prior  approval  of  U.  S.  Army 
Aviation  Materiel  Laboratories,  Tort  Eustis,  Virginia  SH04. 


This  report  presents  the  results  of  a  study  of  the  Influence  of  higher 
order  and  nonlinear  effects  on  the  free  vibration  behavior  of  thin,  cir¬ 
cular  cylindrical  shells.  A  recent  solution  utilizing  the  Karman-Donnell 
strain-displacement  relations  is  examined  and  criticized.  A  new  solution 
is  carried  out  which  removes  the  basis  for  criticism  and  discloses  the 
existence  of  a  nonperiodic  vibration  behavior,  a  phenomenon  heretofore 
unknown.  Further,  solutions  are  obtained  using  the  strain-displacement 
relations  deduced  by  Sanders  and  applied,  in  an  appropriately  modified 
manner,  by  Mayers  and  Rehfield  to  the  shell  postbuckling  problem.  The 
effect  of  the  improved  strain-displacement  relations  in  predicting  the 
vibration  behavior  attendant  to  a  modal  shape  possessing  a  small  number 
of  circumferential  waves  is  assessed,  and  recommendations  are  made  for 
directing  future  effort  on  the  problem. 
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INTRODUCTION 


In  the  past,  the  problem  of  the  instability  of  thin,  circular  cylindrical 
shells  subjected  to  a  variety  of  loadings  has  been  the  object  of  intense 
theoretical  and  experimental  investigation.  A  closely  allied  problem, 
that  of  the  vibrations  of  thin,  circular  cylindrical  shells,  has  also 
come  under  careful  scrutiny,  but  not  to  the  same  extent  nor  with  the 
apparent  success  relative  to  the  instability  problem. 

The  studies  related  to  the  vibrating  cylindrical  shell  problem  have,  for 
the  most  part,  been  contained  within  the  framework  of  linear  theory. 
Notable  contributions  to  the  early  development  of  the  vibrating  cylindri¬ 
cal  shell  problem  are  contained  in  the  work  of  Rayleigh,1  Love,2  FlUgge, J 
and  Arnold  and  Warburton. 4> ^  A  relatively  small  amount  of  work  has  been 
accomplished  concerning  the  vibrational  characteristics  outside  the  realm 
linear  theory  with  the  well-known  Karman-Donnell  equations.  To  assess 
the  merit  of  this  work  and  to  place  the  present  investigation  in  proper 
perspective,  it  is  pertinent  to  dwell  somewhat  on  that  part  of  the  liter¬ 
ature  which  deals  specifically  with  the  free  nonlinear  vibrations  of  thin, 
circular  cylindrical  shells. 


The  first  significant  step  toward  the  -understanding  ana  solution  cf  the 
nonlinear  problem  emanated  from  the  work  of  Eric  Reissner0  in  19pp.  This 
paper,  although  restricted  to  the  linear  theory  of  thin,  shallow  elastic 
shells,  was  the  first  to  show  that  for  vibration  in  which  the  oredominant 
motion  is  radial,  the  equations  of  motion  admit  an  important  simplification. 
Specifically,  for  predominately  radial  motion  of  shallow  shells,  Reissner" 
shows  that  the  midsurface  inertia  terms  (that  is,  the  inertia  terms  due  tc 
motion  in  the  axial  and  circumferential  directions)  can  be  neglected  and 
only  roinor  errors  will  result.  Then,  the  three  coupled  simultaneous  dif¬ 
ferential  equations  of  motion  can  be  simplified  tc  a  system  of  two  simul¬ 
taneous  differential  equations  with  the  dependent  variables  being  the 
well-known  Airy  stress  function  and  the  radial  displacement.  This  latter 
set  of  two  coupled  differential  equations,  independent  of  radial  inertia 
considerations,  becomes  the  well-known  equations  governing  the  static 
instability  problem  of  thin,  circular  cylindrical  shells.  The  basis  for 


1 


this  simplified  set  of  equations  governing  the  phenomena  of  static  in- 

7 

stability  is  contained  in  an  investigation  by  Donnell.  The  work  of 

7 

Donnell,  complementary  to  the  classical  large  deflection  theory  of  flat 
plates,  is  generally  referred  to  today  as  the  Karman-Donnell  shell  theory. 

g 

Later  in  1955,  Reissner  published  a  paper, utilizing  shallow  shell  theory 
dealing  with  both  the  linear  and  nonlinear  nonaxi symmetric  vibrations  of 
thin,  circular  cylindrical  shells.  The  portion  of  the  work  pertaining 
to  the  linear  theory  specializes  the  earlier  work  of  Reissner  to  the 
geometry  peculiar  to  a  thin,  circular  cylindrical  shell.  For  infinitesimal 
deflections,  Reissner  assumes  that  the  radial  displacement  function  is 
representable  as  a  chessboard  pattern.  That  this  assumption  is  valid  is 

strongly  substantiated  by  the  experimental  evidence  of  Arnold  and 

U  5  9  10 

Warburton,  ’  Gottenberg,  and  most  recently  by  Koval  and  Cranch  and 

Koval. In  the  last-noted  reference  of  Koval, ^  no  mention  is  made  of 

comparison  to  linear  theory,  but  a  straightforward  calculation  shows  that 

the  experimental  results  are  almost  identical  to  those  obtained  using  the 

0 

frequency  equation  derived  by  Reissner.  In  the  nonlinear  portion  of 

Q 

Reissner' s  study,  it  is  assumed  that  the  nonlinearity  has  a  more  pro¬ 
nounced  effect  on  the  arbitrary  time  function,  which  modifies  the  choice 
of  deflected  shape,  than  on  the  deflected  shape  itself.  Hence,  the  chess¬ 
board  shape  is  used  also  for  the  nonlinear  vibration  problem.  However, 
even  though  the  chessboard  deflection  pattern  is  a  natural  choice  for 
linear  vibration  problems,  its  selection  in  the  nonlinear  problem  must 
be  more  carefully  assessed.  In  particular,  it  has  been  well  known  in 

shell  postbuckling  problems  (see,  for  example,  the  work  of  von  Karman 
12 

and  Tsien)  that  the  continuity  condition  on  the  circumferential  dis¬ 
placement  must  be  enforced.  This  basic  geometric  constraining  requirement 
must  also  be  employed  in  the  nonlinear  shell  vibration  problem.  However, 
the  use  of  the  chessboard  pattern  leads  to  a  circumferential  displacement 

that  is  not  a  periodic  function  of  the  circumferential  coordinate. 

13 

Evensen  was  the  first  to  point  out  that  this  continuity  requirement 

0 

was  not  satisfied  by  the  choice  of  deflected  shape  assumed  by  Reissner. 
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Another  solution,  similar  to  Reissner's,  was  presented  by  Chu.  Again, 

as  pointed  out  by  Evensen, 1-5  the  continuity  requirement  for  the  circum- 

ferential  displacement  was  not  satisfied.  Evensen  J  further  notes  that 

15 

ir.  a  paper  presented  by  Nowinski,  the  continuity  condition  for  the 
circumferential  displacement  is  taken  into  account;  as  a  consequence, 
however,  a  nonzero  radial  displacement  occurs  at  the  ends  of  the  shell. 

13 

On  the  other  hand,  the  radial  displacement  condition  utilized  by  Evensen, 

8  14 

as  well  as  by  Reissner  and  Chu,  is  that  of  vanishing  radial  displace¬ 

ment  at  the  shell  edges.  It  should  be  noted,  however,  that  this  difference 
in  the  boundary  conditions  would  not  likely  cause  appreciable  difference 
in  the  vibrating  shell  problem  provided  that  the  modal  representation  of 

the  deflected  shape  possessed  a  large  number  of  axial  waves.  This  latter 

13 

point  is  not  mentioned  by  Evensen. 

14  15 

Chu  and  Nowinski  arrived  at  results  which  showed  that  the  nonlinear 
vibrational  characteristics  of  the  shell  exhibited  nonlinearity  of  the 
hardening  type  (that  is,  an  increase  in  the  radial  displacement  magnitude 
is  shown  to  lead  to  an  increase  in  the  frequency  of  radial  vibration) 
which,  in  some  cases,  could  be  strongly  nonlinear.  These  findings  were 
challenged  by  Evensen,1^  who  claimed  that  the  nonlinearity  was  of  the 
softening  type  (that  is,  an  increase  in  the  radial  displacement  magnitude 
is  shown  to  lead  to  a  decrease  in  the  frequency  of  radial  vibration)  and 
was  only  slightly  nonlinear. 

In  the  latter  part  of  1965,  Evensen  and  Fulton1^  presented  a  paper  on  the 
nonlinear  dynamic  response  of  cylindrical  shells  which  is  an  expanded 
version  of  Evensen 's1^  original  investigation.  Evenuen  and  Fulton1^ 
arrive  at  results  which  modify  Evensen' a1 ^  original  conclusion  (that 

is,  that  the  shell  exhibits  either  hard  or  soft  behavior  depending 

upon  the  numerical  values  of  the  wavelength  parameters  attendant  to 

the  assumed  modal  shape).  The  theoretical  results  are  in  agreement 

17 

with  an  experiment  performed  by  Clsur.  in  1965.  However,  the  boundary 
conditions  utilized  in  the  work  of  Evensen  and  Fulton10  reflect  supported, 
rotationally  restrained  ends,  whereas  the  experiment  conducted  by  Clson1^ 
reflects  end  clamping. 


i 


The  objective  of  the  present  investigation  is  to  study  further  the  non¬ 
linear  free  vibrations  of  thin,  circular  cylindrical  shells  on  the  basis 
of  considerations  heretofore  excluded  from  theoretical  work  appearing  in 
the  literature.  It  is  shewn  that  the  energy  level  (potential  plus  kinetic) 
associated  with  a  given  modal  shape  is  an  important  consideration  in  the 
nonlinear  free  vibration  problem.  The  initial  portion  of  the  sti/dy,  based 
upon  the  Karman-Donnell  equations,  reveals  that  the  valid  solution,  cor¬ 
responding  to  a  modal  pattern  possessing  a  large  number  of  axial  waves  and 
described  in  terms  of  free  parameters,  is  the  one  that  reflects  minimum 
energy  considerations,  rather  than  the  one  that  requires  the  vibration 
behavior  to  be  periodic.  This  latter  condition  is  utilized  by  Evensen 
and  Fulton  in  Reference  1 6.  It  is  clearly  demonstrated  that  minimum 
energy  considerations  for  a  free  variation  of  the  modal  pattern  assumed 
in  Reference  1 6  lead  to  a  nonperiodic  motion  for  the  free  nonlinear 
vibrations  of  thin,  circular  cylindrical  shells.  Further,  when  the 
number  of  axial  waves  in  the  modal  pattern  is  small  and  the  influence 
of  boundary  conditions  must  be  accounted  for,  the  adoption  of  a  more 
general  deflected  shape  than  that  employed  in  Reference  16,  but  satis¬ 
fying  the  same  boundary  condition,  again  reveals  the  existence  of  a 
nonperiodic  motion.  In  this  latter  case,  the  nonperiodic  motion  re¬ 
flects  an  energy  level,  for  identical  initial  conditions,  equal  to  that 
associated  with  the  deflected  shape  chosen  in  Reference  16.  The  non¬ 
periodic  motion,  in  all  cases  investigated,  appears  to  be  almost  periodic. 

This  behavior  has  important  implications  concerning  the  interpretation  of 
the  results  obtained  in  experimental  investigations. 

The  present  study  is  not  limited,  as  are  the  previously  referenced  non¬ 
linear  studies,  to  vibrations  where  there  are  a  large  number  of  circum¬ 
ferential  waves  present.  That  there  are  a  large  number  of  circumferential 
waves  is  the  basic  assumption  justifying  the  utilization  of  the  Karman- 
Donnell  equations.  Hence,  to  surmount  the  restriction  on  the  use  of  the 
Karman-Donnell  formulation  of  the  shell  problem  when  the  number  of  cir¬ 
cumferential  waves  is  small,  it  is  necessary  to  employ  a  more  accurate 

« 

set  of  strain-displacement  relations  and,  at  the  same  time,  to  retain 

6  8 

the  midsurface  inertia  terms.  The  observations  by  Reissner  ’  that 

* 
* 

<4 


midsurface  inertia  terms  may  be  neglected  when  the  number  of  circum¬ 
ferential  waves  is  large  can  no  longer  be  invoked. 

The  more  accurate  strain-displacement  relations  employed  are  attributed 

to  Sanders^  and  were  used,  in  a  modified  form,  by  Mayers  and  Rehfield1'* 

in  a  study  of  the  postbuckling  characteristics  of  thin,  circular  cylin- 
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drical  shells.  As  illustrated  by  Mayers  and  Rehfield,  the  postbuckling 
solutions  obtained  on  the  basis  of  the  modified  Sanders  relations  can  be 
utilized  to  establish  the  range  of  validity  of  elastic,  large  deflection 
solutions  (specifically,  those  based  on  the  Karman-Donnell  formulation) 
in  terms  of  shell  radius-to-thickness  ratio.  The  ramifications  of  these 
findings  relative  to  the  nonlinear  shell  vibration  problem  are  pursued. 

Differential  equations  of  motion  and  the  attendant  boundary  terms  are 
developed  on  the  basis  of  a  variational  procedure.  Approximate  solutions 
of  the  equations  are  obtained,  and  the  results  presented  clarify  former 
work  and  yield  new  information  regarding  the  nonlinear  vibrational  char¬ 
acteristics  of  thin,  circular  cylindrical  shells. 
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THEORETICAL  CONSIDERATIONS 


STRAIN-DI3PLACPBTT  RELATIONS 

The  n train -die placement  relations  used  in  this  study  were  deduced  by 

IP 

Sanders.  For  the  thin,  circular  cylindrical  shell,  the  midsurface 
strains  can  be  reduced  to  the  form 
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The  total  strains  are  then  given  by 
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where  ut  ,  v  ,  and  w  (see  Figure  1)  are  the  axial,  tangential,  and 

v  W  w 

radial  cosqjonents  of  displacement,  respectively,  and 
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Substitution  of  Equations  (l)  through  (5)  and  (9)  into  Equations  (6) 
through  (6)  yields  the  strain-displacement  relations 


Equations  (10),  (ll),  and  (12)  are  valid  when  (l)  the  Kirchhoff-Love 
hypothesis  holds  (that  is,  when  a  straight  line  segment  that  is  per¬ 
pendicular  to  the  middle  surface  of  the  undeformed  shell  remains 
perpendicular  to  the  deformed  middle  surface  while  undergoing 
negligible  strain  relative  to  its  original  length),  (2)  the  strains 
and  rotations  of  line  elements  on  the  shell  midsurface  induced  by 
the  deformation  remain  small  in  comparison  to  unity,  although  the 
components  of  the  displacements  are  not  necessarily  small,  and  (3) 
the  angles  of  rotation  about  the  normals  to  the  shell  midsurface  re¬ 
main  negligibly  smaller  than  the  rotations  out  of  the  midsurface 
throughout  the  deformation. 


Although  the  Sanders  Equations  (10),  (11),  and  (12)  appear  as  a  dis¬ 
placement  formulation  in  terms  of  a  tangential  v-displacement  notation, 
they  can  be  written  in  terms  of  a  circumferential  v-displacement  nota¬ 
tion  to  effect  a  convenient  simplification  as  shown  by  Mayers  and 
10 

Rehfield.  '  In  the  alternate  form,  the  equations  become 
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and 
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It  may  be  noted  that  Equations  (lj),  (iL),  and  (15)  reflect  no  quadratic 

terms  in  the  midsurface  displacements  alone.  Also,  for  a  large  number 
of  circumferential  waves,  the  equations  can  be  further  simplified  by 
elimination  of  any  v- displacement  contributions  to  the  curvatures  and 
twist. 


TOTAL  POTENTIAL  ENERGY 

The  total  potential  energy,  in  the  absence  of  body  forces,  consists  of 
the  strain  energy  stored  in  the  body  and  can  be  expressed  as 
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Substitution  of  the  expressions  for  the  strains  from  Equations  (13), 
(14),  and  (15)  into  Equation  (lb)  and  integration  over  the  constant 
thickness  of  the  shell  yield 
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KINETIC  ENERGY 


The  kinetic  energy  is  the  sum  of  the  kinetic  energies  associate!  v:i;5. 
axial,  circumferential,  and  radial  velocities,  respectively;  that  is, 
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Integration  over  the  constant  thickness  yields 
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VARIATIONAL  PRINCIPLE 

The  application  of  Hamilton's  principle  requires  that  the  simultaneous 
first-order  change  in  the  Lagrangian  (L  =  T-U)  ,  integrated  over  a 
specified  time  interval,  with  respect  to  admissiole  variations  in  the 
degrees  of  freedom  characterizing  the  state  of  strain  (namely,  u  ,  v  , 
and  w  )  must  vanish;  that  is, 
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Enforcement  of  this  condition  yields  the  Euler  equations  of  the  variation, 
or  the  differential  equations  of  motion  of  the  shell,  together  with  the 
attendant  boundary  terms. 

EULER  EQUATIONS 

The  Euler  equations  (that  is,  the  displacement  equations  of  motion), 
valid  for  both  periodic  or  nonperiodic  motions,  resulting  from  the 
variation  tvc  given  by 
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The  foregoing  displacement  equations  of  motion  can  be  written  in  a  more 
concise  form  by  employing  Hooke's  law.  Since  Hooke's  law  states  that 
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then,  though  utilization  of  Equation  (2k),  the  equations  of  motion  can 
be  written  as 
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The  equations  can  be  written  in  a  still  more  compact  form  by  the  intro¬ 
duction  cf  shear  ar d  bending  moment  resultants  defined  as  fellows  (see 
Figure  2' : 
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Then  use  of  the  definitions  given  by  Equation  (28)  allows  the  equations 


of  motion  to  be  written  as 
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or  in  the  linearized  form  as 
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The  linearized  equations  of  notion,  reduced  to  the  static  case,  are  not 
identical  to  the  linear  set  presented  by  Sanders  in  Reference  20.  The 
only  difference  between  the  two  sets  of  linear  equations  ic  a  consequence 
of  the  neglect,  in  the  present  development,  of  effects  associated  with 
rotations  about  the  normal  to  the  shell  midsurface. 


BOUNDARY  CONDITIONS 

The  boundary  conditions  which  result  from  application  cf  Hamilton's 
principle  (Equation  (20))  are  given,  for  both  periodic  and  nonperiodic 
motions,  in  general  form  as 

hcr^  =  0  or  5u  =  0  at  x  =  C  ,l  (35) 
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REDUCTION  TO  CLASSICAL  EQUATIONS 


A  simplified  set  of  Euler  equations  can  be  obtained  if  the  Karman-Donnell 
strain-displacement  relations  are  used  in  place  of  the  more  complicated 
set  of  strain-displacement  relations  given  by  Equations  (13),  (1*0,  and 
(15).  The  Karman-Donnell  strain-displacement  relations  are  given  by 
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Substituting  Equations  (39),  (40),  and  (41)  into  the  Lagrangian  and 
carrying  out  the  variations  in  accordance  with  Hamilton's  principle  - 
leads  to  the  following  simplified  set  of  stress-displacement  equations 
of  motion  for  the  axial,  circumferential,  and  radial  directions, 
respectively. 
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Now,  if  the  notion  ia  predominantly  radial  (as  when  there  are  a  large 
number  of  circumferential  waves  in  the  modal  pattern),  then  Equations 
(42),  (1*5),  and  (44)  admit  still  another  and  quite  drastic  simplifica¬ 
tion,  that  cf  neglecting  the  aidaurface  inertia  terms.  The  Justification 
for  this  simplification  was  first  r.ctcd  ny  ReissndT.  Thus,  with  the 
midsurface  inertia  terms  omitted,  the  equilibrium  eq-atiqns  become 
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Equations  (45)  ar.d  (46)  car.  now  be  satisfied  identically  by  the  intro¬ 
duction  of  ar.  Airy  stress  function  F  defined  such  that 
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Then,  use  of  Equations  (4-5),  (46),  and  (48)  permits  the  lateral 
equilibrium  equation  to  be  written  as 
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A  second  equation  involving  the  stress  function  and  the 
ment  can  be  obtained  by  combining  Equations  (24),  (39), 


(46),  and  (48)  to  yield 
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This  compatibility  equation  was  developed  first  by  Donnell.  Equations 
(49)  and  '50),  with  the  radial  inertia  term  omitted,  are  one  of  the 
equation  sets  commonly  referred  to  as  the  Karman-Donnell  equations; 
these  equations  govern  the  large-deflection  behavior  of  thin  shells  in 
the  presence  of  midsurface  forces. 
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METHOD  OF  SOLUTION 


The  differential  equations  obtained  in  the  previous  section  which  re¬ 
flect  the  utilization  of  strain-displacement  relations  equivalent  to 
those  of  Sanders  are  three  highly  nonlinear  coupled  equations  in  the 
three  dependent  variables  u  ,  v  ,  and  w  .  On  the  other  hand, 
Karman-Donnell  strain-displacement  relations,  with  the  midsurface 
inertia  terms  neglected,  admit  the  reduction  to  two  simultaneous 
differential  equations  involving  only  the  stress  function  and  the 
radial  displacement.  This  reduced  set  of  two  simultaneous  equations 
is  still  highly  nonlinear  and  coupled  but  is  much  more  manageable  than 
the  set  of  three  simultaneous  differential  equations.  The  solutions 
reflecting  the  use  of  both  sets  of  strain-displacement  relations  are 
obtained  by  employing  the  direct  method  of  the  calculus  of  variations 
(that  is,  the  Rayleigh-Ritz  method).  The  solutions  corresponding  to 
the  strain-displacement  relations  which  result  in  the  Karman-Donnell 
equations  ((49)  and  (50))  are  discussed  first  for  the  case  when  a  large 
number  of  circumferential  waves  are  present.  The  solutions  corresponding 
to  the  Sanders  strain-displacement  relations  (Equations  (21),  (22),  and 
(23))  are  discussed  subsequently  for  the  case  when  a  small  number  of 
circumferential  waves  occurs. 

SOLUTIONS  FOR  A  LARGE  HUMBER  OF  CIRCUMFERENTIAL  WAVES 

In  all  of  the  existing  analytical  studies  of  the  nonlinear  vibrations  of 
thin,  circular  cylindrical  shells,  except  for  the  work  reported  in  Ref¬ 
erences  13  and  l6,  it  is  assumed  that  (l)  the  radial  displacement  is 
representable  by  a  chessboard  pattern,  (2)  the  shell  possesses  many 
circumferential  waves  (that  is,  the  shallow-shell  behavior),  and  (3)  the 
midsurface  inertia  terms  can  be  neglected.  If  indeed  there  are  a  large 
number  of  circumferential  waves  present,  then  assumption  (3)  seems  in¬ 
tuitively  correct.  Assumption  (1),  however,  seems  intuitively  incorrect. 
As  noted  by  References  13  and  16,  utilization  of  the  chessboard  pattern 
for  the  radial  displacement  leads  to  a  violation  of  the  periodicity  con¬ 
dition  for  the  circumferential  displacement.  Even  if  the  periodicity 
condition  relative  to  a  chessboard  pattern  for  the  radial  displacement 
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could  be  satisfied  by  selecting  a  special  homogeneous  solution  for  the 
stress  function  (Reference  13  alludes  to  having  accomplished  this,  but 
no  details  are  presented),  it  still  seems  intuitively  incorrect  that  the 
radial  displacement  pattern  for  a  large  number  of  circumferential  waves 
would  be  am  exact  chessboard  pattern.  Some  consideration  must  be  given 
to  the  degree  of  midsurface  stretching  and,  consequently,  to  the  increased 
level  of  energy  that  would  be  necessary  to  produce  a  pure  chessboard  pat¬ 
tern.  It  would  appear  that,  for  finite  radial  displacements,  a  displace¬ 
ment  pattern  which  allows  the  shell  to  move  inward  would  be  more  likely 
to  occur  since  this  type  of  displacement  pattern  should  involve  less 
energy  and,  hence,  lead  to  the  equilibrium  position  that  the  shell  would 
seek.  This  intuitive  reasoning,  which  is  founded  on  a  knowledge  of  shell 
behavior  in  the  axial  compression  postbuckling  problem,  is  shown  in  retro¬ 
spect  to  be  justified. 

Fixed  Parameter  Solution  (n  Large) 

For  a  uniform  homogeneous  shell,  the  study  reported  in  Reference  1 6  is 
based  on  the  assumption  that  the  radial  displacement  which  satisfies  the 
geometrical  constraints  of  the  problem  (that  is,  freely  supported  ends) 
is  representable,  save  for  a  shift  of  coordinate  system  origin,  in  the 
following  form: 

(w  \  nx  ity  T)  2  2nx  rj  2 

-  )  =  yn  (t)  cos  —  cos  —  +  -  y.  (t)  cos  - +  -  y.  (t)  (51) 

\  \  8  8 

where  y.  (t)  is  an  arbitrary  function  of  time.  The  above  choice  of  de¬ 
flected  shape,  while  not  being  perfectly  general,  is  a  discerning  one 
in  that,  in  addition  to  satisfying  the  geometric  boundary  conditions, 
it  also  leads  to  satisfaction  of  the  continuity  constraint  on  the  cir¬ 
cumferential  displacement.  Further,  by  restricting  the  time  dependent 
coefficients  to  terms  proportional  to  y^t)  and  y  (t)  ,  the  resultant 

motion  is  forced  to  be  periodic,  although  not  harmonic. 

In  Reference  16,  the  Galerkin  procedure  was  used  to  obtain  the  governing 
differential  equations;  subsequent  solution  of  the  equations  was  carried 
out  by  means  of  a  perturbation  technique.  Since  the  Galerkin  integral  ~ 
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Approach  ll.-.closes  no  information  relative  to  energy  level,  the  present 
investigation  reformulates  and  solves  the  problem  by  using  the  vari¬ 
ational  approach:  the  radial  displacement  given  by  Equation  (51),  to¬ 
gether  with  the  stresses  obtained  from  Equations  (48)  and  (50),  is 
substituted  into  the  Lagrangian  (see  Appendix  I  for  development).  The 
enforcement  of  the  vanishing  of  the  first  ve  s  t  on  of  the  Lagrangian 
with  respect  to  y  (t)  leads  to  a  single  sl  ond-order  nonlinear  dif¬ 
ferential  equation  with  y^(f)  as  dependent  variable.  After  re¬ 
duction  of  the  single  second-order  differential  equation  to  two  coupled 
first-order  differential  equations,  the  solution  is  obtained  by  employing 
a  standard  Kutta-Merson  numerical  method. 

The  results  obtained  verify  that  the  curves  in  Reference  16  representing 
the  nonlinear  free,  but  periodic,  vibrational  characteristics  are  indeed 
correct.  These  results  are  shown  in  Figures  3  and  4. 

However,  the  occurrence  of  the  multivalued  conditions  of  (1)  different 
amplitudes  for  a  given  wave  pattern  at  a  fixed  value  of  the  frequency 
and  (2)  different  wave  patterns  for  the  same  amplitude  at  a  given  fre¬ 
quency  requires  the  approach  of  the  present  analysis,  which  establishes 
the  lowest  level  of  energy  associated  with  a  particular  condition.  The 
discussion  of  this  point  is  expanded  in  the  section  on  "Results  and 
Discussion".  It  is  also  observed  that  the  periodic  solution,  which  is 
based  upon  a  vanishing  radial  deflection  at  the  cylinder  ends,  implies 
that  the  vibrational  behavior  is  influenced  by  the  boundary  conditions 
in  the  presence  of  a  modal  pattern  containing  a  small  number  of  axial 
waves.  Although  the  case  of  a  small  number  of  axial  waves  is  not  the 
general  problem  being  investigated  herein,  it  will  be  shewn  in  the  next 
section  that  a  more  general  radial  deflection  function  (one  which  results 
in  the  vanishing  of  the  radial  deflection  at  the  shell  ends  and  leads  tc 
satisfaction  of  the  v-displacement  periodicity  condition)  does  not  reflect 
necessarily  a  periodic  motion  based  upon  the  minimum  energy  criterion. 

For  long  shells  (large  number  of  axial  waves),  the  nonperiodic  vibrational 
behavior  is  shown,  for  all  cases,  to  correspond  to  the  minimum  energy 
solution  when  the  parameters  of  the  radial  deflection  function  of  Reference 
16  are  free  rather  than  fixed. 
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Free  Parameter  Solution  (n  Large) 

If  the  shell  modal  pattern  reflects  several  or  more  waves  in  the  axial 
direction,  then  it  is  reasonable  to  assume  that  the  boundary  conditions 
will  not  significantly  influence  the  frequency  of  vibration.  Then  the 
deflected  shape  used  in  the  fixed  parameter  solution  of  the  previous 
section  can  be  made  more  general  and  representable  as 

(wv  nx  ny  2nx 

-)  =  y-,(t)  cos  —  cos  —  +  y,(t)  cos  -  +  yslt)  (52) 

h  /  IX'  X 

x  y  x 

where  y^(t)  ,  y.(t)  ,  and  y^t)  are  arbitrary  functions  of  time.  The 
substitution  of  the  radial  displacement  given  by  Equation  (52)  and  its 
associated  stress  function,  obtained  from  Equation  t.50),  into  the  con¬ 
tinuity  condition  for  the  circumferential  displacement  yields  the  re¬ 
quirement 


y5(t)  =  Jy^t) 


(53) 


Hence,  the  deflected  shape  is  now  given  by 

ny 


nx 

(t)  cos  —  cos 
X 


2nx  T]  2 

+  y  ,(t)  cos  -  +  -  y,  (t) 

\  j  X  8  1 

y  x 


(54) 


The  substitution  of  the  radial  displacement  given  by  Equation  (54),  to¬ 
gether  with  the  stresses  obtained  from  Equations  (48)  and  (50),  into  the 
Lagrangian  (see  Appendix  I  for  development),  and  enforcement  of  the 
vanishing  of  the  first  variation  of  the  Lagrangian  with  respect  to 
y^(t)  and  y^(i)  lead  to  two  second-order  nonlinear  differential 
equations  with  y^(x)  and  y^(t)  as  the  independent  variables.  After 
a  change  of  the  independent  variable  and  subsequent  reduction  of  the 
two  second-order  differential  equations  to  four  first-order  differential 
equations,  the  solution  is  obtained  by  employing  a  standard  Kutta-Merson 
numerical  method. 


A  comparison  of  the  maximum  radial  displacement  for  the  fixed  parameter 
solution  possessing  the  same  modal  shape  and  the  same  initial  conditions 
is  shown  on  Figure  5.  The  free  parameter  solution,  reflecting  only  a 
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change  in  the  initial  value  of  y,  , 
solution  is  shown  in  Figure  6. 


leading  tc  the  minimum  energy 


A  comparison  of  the  nondimensional  strain  energy-  parameter  for  the  cases 
where  the  initial  velocities  are  zero  is  given  in  Figures  7,  ’’ ,  and  9 
for  different  combinations  of  the  modal  shape  and  initial  displacements. 


When  the  shell  modal  pattern  reflects  a  small  number  of  axial  waves,  the 
choice  of  deflected  shape,  which  enforces  the  vanishing  radial  deflection 
condition  at  the  shell  ends,,  can  be  made  more  general  than  that  given  in 
Reference  16  and  is  representable  as 


2nx  t]  „ 

-  y,(t)  cos  -  +  -y,  (t)  (55) 

3  \  e 

The  radial  deflection  (Equation  (55))  leads  tc  identical  satisfaction 
of  the  v-displacement  periodicity  condition.  The  resulting  solution, 
obtained  in  a  manner  analogous  to  the  previous  free  parameter  solution, 
yields  a  nonperiodic  vibrational  behavior.  A  comparison  of  the  maximum 
radial  displacements,  reflecting  the  different  assumed  deflection  func¬ 
tions,  is  given  on  Figure  5. 


SOLUTION  FOR  A  SMALL  NUMBER  OF  CIRCUMFERENTIAL  WAVES 

For  a  uniform  homogeneous  shell,  the  strain  energy  and  the  kinetic  energy 
are  given  by  Equations  (.17)  and  (10),  respectively.  Functions  representing 
the  degrees  of  freedom  u  ,  v  ,  and  w  ,  which  satisfy  the  continuity  con¬ 
dition  on  the  circumferential  displacement  but  are  arbitrary  with  time, 
may  be  taken  in  the  form 


(•)■ 


nx  jty 
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x  y 
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Jtx  jty  2itx  2jty 

=  •  y  (t)  cos  —  sin  —  +  y  (t)  cos  -  sin  -  +  y,  .(t) 
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where  y1(t)  ,  y^(t)  ,  y^(t)  ,  y^(t)  ,  y^(t)  ,  yu(t)  ,  and  y^Ct)  are 
arbitrary  functions  of  time.  Now,  the  substitution  of  Equation  (56)  in¬ 
to  Equations  (17)  and  (18)  and  the  enforcement  of  the  simultaneous  vanish¬ 
ing  of  the  first  variation  of  the  Lagrangian  (L  =  T-U)  with  respect  to 
y1(f)  ,  y3(f)  ,  y5(t)  ,  y  (t)  ,  y9(T)  ,  y^T)  ,  and  y^t)  le«i  to  a 
set  of  seven  second-order  nonlinear  coupled  differential  equations. 

After  a  change  of  the  independent  variable  and  subsequent  reduction  of 
the  7  second-order  differential  equations  to  lU  first-order  differential 
equations  (the  details  of  which  are  given  in  Appendix  III),  the  solution 
is  obtained  by  employing  a  standard  Kutta-Merson  numerical  method. 


The  degree  and  character  of  the  nonlinearity,  reflecting  the  employment 
of  the  modified  Sanders  strain-displacement  relations,  are  shown  in 
Figures  10  and  11  for  the  case  of  finite  initial  displacements  but  zero 
initial  velocities. 
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RESULTS  AMD  DISCUSSION 


In  the  free  vibration  of  a  physical  system  represented  or  idealized  by 
a  set  of  linearized  equations,  the  magnitude  of  the  resulting  displace¬ 
ments  and  the  energy  associated  with  these  displacements  is  indeterminate. 
However,  if  the  equations  representing  the  physical,  system  are  not  lin¬ 
earized,  the  displacements  and  the  attendant  energy  then  possess  a  unique 
determinate  value  for  any  given  set  of  initial  conditions. 

It  is  known  that  in  the  solution  of  the  postbuckling  problem  of  a  thin, 
circular  cylindrical  shell,  the  determination  of  the  minimum  energy  is  of 
primary  importance.  One  of  the  solutions  given  in  Reference  12  predicts 
that  the  shell  will  be  in  tension  in  the  postbuckled  state  under  the  ac¬ 
tion  of  an  applied  compressive  end  load.  This  predicted  solution  is 
obviously  incorrect,  as  is  pointed  out  in  Reference  21.  The  solution  is 
rejected  because  it  does  not,  in  general,  represent  the  minimum  energy 
for  the  particular  choice  of  radial  deflection  function  selected.  Hence, 
since  energy  considerations  are  important  in  the  understanding  and  compre¬ 
hension  of  the  nonlinear  shell  postbuckling  problem,  it  would  appear  that 
these  same  considerations  would  prove  to  be  vitally  important  to  any  solu¬ 
tions  of  the  nonlinear  shell  vibration  problem. 

That  the  energy  level  for  a  given  modal  shape  might  be  an  important  con¬ 
sideration  has  not  been  brought  out  or  explored  to  date.  The  most  advanced 
solution  of  the  nonlinear  shell  vibration  problem  (Reference  16)  makes  no 
mention  of  energy  considerations  relative  to  the  assumed  modal  shape.  It 
is  stated  in  Reference  lU  that,  for  a  large  number  of  circumferential  waves, 
the  square  wave  pattern  (that  is,  the  chessboard  pattern)  has  the  least 
strain  energy.  However,  this  statement  is  not  correct,  and  this  point  is 
discussed  later  in  this  section. 
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Still  another  important  facet  in  the  understanding  of  the  nonlinear  free 
vibration  problem  of  a  thin,  circular  shell  involves  an  analogy  to  the 
static  instability  problem.  In  the  nonlinear  shell  postbuckling  problem 
(see.  for  example,  Reference  22),  the  simultaneous  minimization  is  made 
with  respect  to  the  free  parameters  in  the  assumed  deflected  shape  and 
also  with  respect  tc  the  modal  wave  parameters  u  and  q  .  As  such, 
the  number  of  axial  and  i,ircumferential  waves  is  continuously  varying  as 
the  minimization  is  accomplished.  However,  the  shell  instability  problem 
is  treated  by  means  of  a  static  analysis,  and  a  slightly  different  approach 
is  necessary  in  the  free  vibration  problem.  In  the  linear  and  nonlinear 
free  vibration  problem,  the  basic  assumption  is  introduced  that  the  wave 
parameters  u  and  ^  are  net  variables  but  are  constant.  This  assump¬ 
tion  essentially  implies  that,  if  an  experiment  were  being  conducted,  and 
if  the  modal  shape  and  modal  lines  (w  =  0)  could  be  detected,  then  this 
modal  shape  and  its  modal  lines  would  not  vary  with  time.  The  implication 
of  this  assumption  relative  to  the  assumed  modal  shapes  chosen  in  the  pre¬ 
vious  section  is  discussed  subsequently. 

The  deflected  shape  chosen  in  Reference  l6  (that  is,  the  fixed  parameter 
solution  reflecting  a  large  number  of  circumferential  waves)  is  a  judicious 
choice.  It  has  been  shown  in  References  h,  5,  9,  10,  and  11  that  the  chess¬ 
board  pattern  is  indeed  the  pattern  which  the  shell  seeks  when  the  displace¬ 
ments  involved  are  infinitesimal.  Hence,  it  would  seem  logical  to  assume 
that  a  slightly  modified  chessboard  pattern  might  emerge  as  the  modal  pat¬ 
tern  when  the  displacements  are  increased  from  infinitesimal  to  finite 
values.  In  the  analogous  shell  buckling  problem  (see,  for  example  Ref¬ 
erences  12  or  22),  three  of  the  most  important  terms  in  the  representation 
of  the  radial  deflection  are  those  selected  in  Reference  16.  In  addition, 
the  deflected  shape  selected  in  Reference  16  reflects  zero  radial  displace¬ 
ment  at  the  ends  of  the  shell  and  leads  to  a  periodic  solution  in  time  when 
the  continuity  condition  on  the  v-displacement  is  enforced. 
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Figures  3  and  4  depict  the  frequency  behavior  for  the  resultant  periodic 
motion  as  a  function  of  the  amplitude.  The  nonlinearity  of  the  resulting 
motion  and  the  degree  of  hardness  or  softness  of  the  system  are  highly 
dependent  upon  the  value  of  the  wave  parameter  q  ;  for  the  same  initial 
displacement  (assuming  an  initial  velocity  of  zero),  the  nonlinearity 
becomes  more  pronounced  as  the  value  of  t)  increases.  Hence,  for  a 
particular  shell  geometry  (fixed  (h/R)  )  and  a  fixed  set  of  initial 
conditions,  the  degree  of  nonlinearity  is  increased  as  the  number  of 
circumferential  waves  becomes  larger. 

It  can  be  further  noted  from  Figures  3  and  4  that  the  character  of  the 
nonlinearity  is  dependent  upon  the  value  of  the  full  wave  aspect  ratio 
p  .  In  general,  the  results  show  that  the  vibrations  usually  exhibit 
a  softening  character  when  p  <  1  and  a  hardening  character  when  p  >  1  . 
In  the  present  analysis,  the  strain  energies  associated  with  the  various 
choices  of  initial  conditions  and  wave  parameters  have  been  determined. 
Several  values  of  the  energy  appear  on  Figures  3  and  4.  In  particular, 
it  is  noted  that  at  a  certain  set  of  fixed  initial  conditions  (an  initial 
displacement  of  unity  and  air.  initial  velocity  of  zero)  the  energy  increases 
considerably  as  the  value  of  the  full  wave  aspect  ratio  p  increases  (see 
Figure  4).  This  is  more  clearly  shown  in  Figure  12.  In  a  like  manner,  the 
increase  in  energy  associated  with  an  increase  in  the  wave  parameter  q  is 
shown  on  Figure  13. 

It  is  noted  that  at  the  crossover  points  on  Figure  3,  the  energies  are  not 
the  same.  For  instance,  at  the  point  corresponding  to  the  coordinates 
y1  =  1.91  and 

modal  shape  is  185#  higher  than  the  energy  associated  with  the  p  =  0.707 
modal  shape.  This  indicates  that  for  the  same  frequency  and  displacement, 
the  modal  shape  is  dependent  upon  the  energy.  Hence,  in  am  experiment 
where  a  large  frequency  spectrum  is  examined,  certain  modes  may  never 
appear  if  the  vibration  excitor  camnot  inject  sufficient  energy  to  the 
shell.  This  is  particularly  true  for  the  higher  levels  of  energy  which 
occur  for  p  >  1  . 
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Another  point  not  mentioned  in  Reference  16,  but  worthy  of  mention  here, 
can  be  noted  in  Figure  14,  which  is  a  partial  plot  of  the  modal  shape 
for  the  listed  choice  of  deflected  shape.  As  can  be  seen  in  Figure  14, 
the  aspect  ratio  is  not  constant  but  is  variable.  However,  if  the  whole 
wavelength  of  modal  shape  is  considered  and  the  aspect  ratio  of  the  wave¬ 
form  is  taken  as  the  average  ever  this  whole  wavelength,  then  it  can  be 
considered  as  a  constant. 

Further,  it  should  be  pointed  out  that  the  solution  obtained  in  Reference 
lo  is  presumably  valid  for  the  range  of  values  of  the  wave  parameters 
presented.  As  such,  the  solution  should  be  valid  when  the  number  of 
axial  waves  is  large  or  small.  However,  it  appears  that  while  the 
solution  presented  in  Reference  1 6  is  the  most  advanced  solution  to 
date,  new  insight  into  the  problem  is  gained  through  energy  considera¬ 
tions.  The  result  of  freeing  all  parameters  in  the  assumed  radial  de¬ 
flection  fractions,  compatible  with  assuring  the  continuity  condition 
on  the  v-displacement,  leads  to  nonperiodic  motions  when  the  number  of 
axial  waves  is  either  large  or  small.  Also,  although  the  motion  is  non¬ 
periodic,  strictly  speaking,  in  certain  cases  it  appears  similar  to  the 
periodic  fixed-parameter  solution  (see  Figures  5  and  6).  Hence,  too 
casual  an  observation  of  experimental  results  might  lead  to  the  erroneous 
conclusion  that  the  motion  is  periodic. 

*  -  -  P  * 

When  the  number  of  axial  waves  is  small,  the  energy  levels  reflecting  the 
deflected  shapes  (Equations  (51)  or  (55))  are  identical  for  the  present 
choice  of  initial  conditions,  but  the  more  general  choice  of  deflected 
shape  (Equation  (55))  leads  to  a  nonperiodic  motion  while  the  solution 
of  Reference  1 6  (that  is,  Equation  (51))  leads  to  a  periodic  motion  (see 
Figure  5).  Although  the  case  of  a  small  number  of  axial  waves  is  not  the 
case  of  interest  for  the  present  investigation,  it  is  highly  probable  that 
introducing  a  still  more  general  radial  deflection  function  for  the  "freely 
supported"  case  will  yield  a  nonperiodic  solution  which  possesses  a  lower 
energy  level  than  that  associated  with  the  radial  deflection  functions 
given  by  Equations  (51)  and  (55).  This  conjecture  is  based  on  the  known 
tendency  for  thin  cylinders  deformed  into  many  waves  to  depart  more  and 
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more  from  ohe  chessboard  modal  pattern  as  larger  and  larger  amplitudes 
of  the  radial  deflection  occur. 

When  the  number  of  axial  waves  is  large,  the  nonperiodic  solution,  which 
reflects  the  radial  deflection  given  by  Equation  (51*-),  leads  to  a  lower 
energy  than  the  periodic  solution  (Equation  (51)).  The  lower  energy, 
associated  with  the  nonperiodic  solution,  can  be  seen  on  Figures  7,  8, 
and  9.  On  Figure  7,  the  minimum  energy  occurs  almost  identically  at 
the  value  of  y^  which  Reference  16  employs.  That  is,  however,  only 
chance  occurrence.  In  Figure  8,  the  true  minimum  energy  solution  re¬ 
flects  a  significant  difference  in  the  value  of  y^  from  that  which 
was  utilized  in  Reference  16.  Further,  it  is  interesting  to  note  that 
a  lower  value  of  the  y^  displacement  shifts  the  y^  value  for  the 
minimum  energy  from  a  higher  to  a  lower  value.  From  this  observation 
and  a  perusal  of  the  strain  energy  expression  for  the  free  parameter 
solution  (see  Equation  (68)),  an  important  physical  consequence  can  be 
deduced.  As  the  value  of  y^  becomes  smaller  and  smaller  (that  is,  as 
it  approaches  infinitesimal  values),  the  minimum  point  shifts  farther 
and  farther  to  the  left;  this  leads  to  the  conclusion  that  for  infinitesimal 
values  of  y^  ,  the  minimum  energy  is  the  solution  corresponding  to 

y,  =  0  .  This  observation,  coupled  with  the  fact  that  y1  <  1  implies 

J  2  x  ■ 

that  y^  «  1  ,  leads  to  the  conclusion  that  in  approaching  the  limit 

of  infinitesimal  displacements,  the  least  energy  solution  reflects  the 

chessboard  pattern.  That  this  is  true  can  also  be  deduced  from  Figure 

This  figure  shows  that  for  a  small  value  of  displacement  (y^  =  0.01)  , 

the  minimum  energy  solution  exists  for  y^  =  0  ;  hence,  further  credence 

is  given  to  the  deduction  discussed  above.  This  deduction  is  important. 

It  shows  that  only  when  the  displacements  are  very  small  is  the  minimum 

energy  solution  representable  as  a  chessboard  pattern;  it  also  shows  that 

for  finite  displacements  of  the  order  of  magnitude  of  the  shell  thickness, 

the  minimum  energy  is  related  to  a  slightly  modified  chessboard  pattern. 

This  shows  why  in  the  linear,  small  displacement  experimental  studies 

(References  U,  5,  9,  10,  and  11),  the  predominant  evidence  indicated  that 

the  chessboard  pattern  was  the  actual  pattern  developed  by  the  vibrating 

shell. 
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Still  another  and  important  consideration  deals  with  the  validity  of 
the  solution  concerning  the  rotation  limits  allowable  within  the  scope 
of  Karman-D 'nnell  theory.  'This  consideration  was  used  in  Reference  19 
in  discussing  limitations  on  the  validity  of  several  well-known  solu¬ 
tions  of  the  postbuckled  shell  problem.  The  maximum  rotations  must  be 
small  in  comparison  to  ’unity  (5  percent  is  considered  herein  to  be  the 
limiting  value).  Then,  Figures  15,  16,  17,  and  18  show  conservative 
estimates  of  the  limitations  inherent  in  the  solution.  The  severe 
limitations  on  the  p  =  2  case  (Figure  18)  should  be  noted.  For  a 
y^  =  4  ,  the  R/h  of  the  shell  must  be  greater  than  or  equal  to  1-500 
to  be  within  the  range  of  allowable  rotations.  Thus,  for  shells  which 
are  in  the  range  of  practical  interest  (100  <  —  <  3000)  ,  the  validity 

of  the  present  solution  would  have  to  be  restricted  (for  the  u  •=  2 
case)  to  smaller  amplitudes.  It  may  be  noted  further  that  the  restric¬ 
tions  on  the  validity  of  the  solutions  are  less  stringent  for  other 
values  of  n  (see  Figures  15,  16,  and  17).  Hence,  unless  the  displace¬ 
ments  experienced  by  the  shell  are  relatively  large,  the  limitations 
concerning  the  allowable  rotation  do  not  appear  to  be  critical  for  the 
free  vibration  problem.  However,  these  limitations  could  become  con¬ 
siderably  more  important,  for  example,  in  a  dynamic  buckling  or  a  forced 
vibration  study. 

For  vibration  with  low  values  of  n  ,  the  Sanders  strain-displacement 
relations  have  been  used  and  the  midsurface  inertia  terms  retained.  With 
the  retention  of  the  midsurface  inertia  terms,  the  occurrence  of  three 
frequencies  for  each  given  modal  configuration  is  possible  (see,  for 
example,  References  4  and  5).  However,  each  of  the  frequencies  will 
have  different  amplitude  ratios;  of  the  three  possible  values  of  the 
frequency,  the  important  one  is  that  associated  with  the  primarily 
radial  motion  of  the  shell  walls.  This  value  is  several  orders  of 
magnitude  lower  than  either  of  the  other  two,  which  are  associated 
primarily  with  axial  and  circumferential  motions,  respectively.  The 
deflected  shapes  and  the  attendant  initial  conditions  assumed  in  this 
investigation  again  reflect  the  interest  in  only  the  motion  which  is 
predominantly  radial. 
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Due  to  the  complexity  of  this  nonlinear  shell  vibration  problem,  the 
free  vibration  solutions  presented  herein  are  based  only  upon  the  basic 
chessboard  modal  pattern.  For  small  amplitudes,  the  cness board  pattern 
continues  to  provide  a  periodic  motion,  which  is  not  necessarily  the 
true  situation.  Thus,  where  quantitative  rather  than  qualitative  behavior 
is  of  prime  importance,  the  chessboard  pattern  must  be  modified  to  permit 
the  possibility  of  the  occurrence  of  a  nonperiodic  vibration  behavior. 

In  the  modal  pattern  associated  with  n  =  2  ,  the  nonlinearity  of  the 

resultant  motion,  as  in  the  case  when  the  number  of  circumferential- 
waves  is  large,  is  highly  dependent  upon  the  value  of  the  wave  parameter 
q  .  For  a  fixed  number  of  circumferential  waves  (here  n  =  2  ) ,  the 

low  values  of  q  reflect  the  behavior  of  a  thin  shell,  whereas  the  larger 
values  of  q  reflect  the  behavior  of  a  thicker  shell.  Representative 
curves  depicting  this  behavior  are  shown  in  Figures  10  and  11.  In 
Figures  10  and  11,  it  is  shown  that  no  detectable  nonlinearity  is 
observed  until  the  (R/h)  ratio  is  less  than  1000  .  For  gradually 
increasing  thickness,  the  nonlinearity  becomes  pronounced.  In  Figure  10 
for  a  full  wave  aspect  ratio  of  unity,  the  ratio  of  nonlinear  to  linear 
frequency  for  (h/R)  =  1/100  and  (w/h)  =  16  reaches  1.165  or  16.5$ 
beyond  that  which  linear  theory  would  predict.  In  Figure  11  for  a  full 
wave  aspect  ratio  of  1/2  ,  the  ratio  of  nonlinear  to  linear  frequency, 
again  for  (h/R)  =  1/100  and  (w/h)  =  16  ,  reaches  1.25  or  25$  over 

that  which  linear  theory  predicts. 


For  the  large  values  of  deflection  necessary  to  produce  a  pronounced 
nonlinearity  when  the  number  of  circumferential  waves  is  low,  it  is 
imperative  to  ascertain  at  what  point  the  limitations  on  the  rotations 
are  exceeded.  Representative  values  for  q  =  0.0k  and  u  =  1.0  are 
shown  on  Figures  19  and  20.  While  these  limitations  are  intended  to 
be  only  a  qualitative,  conservative  estimate,  it  is  apparent  that  the 

In 

,  the  data  of  Figure  20  predict 


2  2 
(cp^)  limitation  is  more  critical  than  the  (cp^)  limitation. 


particular,  for  q  =  O.Ok  and  pi  =  1.0 
2 

that  (cp^)  will  exceed  its  allowable  limit  for  an  (R/h)  =  100  some¬ 
where  between  the  point  where  the  initial  y.  displacement  is  between 
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10  and  12.  For  y^  =  12  ,  the  character  of  the  nonlinearity  is  given 

in  Figure  10.  Again,  the  consideration  of  the  size-of-rotations  limita¬ 
tion  on  the  validity  of  the  theory  does  not  appear  to  be  critical  for 
the  free  vibration  problem.  However,  large  amplitude  behavior,  in  con¬ 
nection  with  dynamic  buckling  or  forced  vibration  analyses,  could  lead 
to  restrictions  on  the  validity  of  solutions  for  specific  shell  geometries. 


CONCLUDING  REMARKS 


The  governing  differential  equations  and  the  attendant  boundary  conditions 
for  the  nonlinear  free  vibration  of  thin,  circular  cylindrical  shells  have 
been  developed  by  variational  procedures  by  the  use  of  both  the  Karman- 
Donnell  and  a  set  of  appropriately  modified  Sanders  strain-displacement 
relations.  Approximate  solutions,  obtained  by  the  direct  variational 
approach,  are  utilized  to  re-examine  and  provide  new  interpretation  to 
established  approximate  solutions  for  shells  possessing  a  large  number 
of  circumferential  waves  (that  is,  when  the  Karman -Donnell  formulation 
is  valid).  Approximate  solutions,  which  reflect  the  use  of  the  modified 
Sanders  formulation  for  shells  possessing  a  small  number  of  circumferential 
waves,  are  obtained  and  are  believed  to  be  the  first  application  of  this 
higher  order  shell  theory  to  a  nonlinear  free  vibration  problem. 

The  solutions,  which  reflect  the  utilization  of  the  Karman -Donnell  for¬ 
mulation  for  a  representative  assumed  modal  shape  possessing  a  large  or 
small  number  of  axial  waves,  reveal  the  existence  of  nonperiodic  motions. 
When  it  is  reasonable  to  assume  that  the  boundary  conditions  at  the  shell 
ends  will  not  significantly  affect  the  frequency  of  vibration  (that  is, 
when  the  modal  shape  possesses  a  large  number  of  axial  waves),  the  non¬ 
periodic  solution  reflects  a  lower  energy  level  than  the  established 
periodic  solution  which  always  implies  the  boundary  condition  of  "freely 
supported"  ends.  When  the  boundary  conditions  have  a  more  pronounced 
effect  on  the  vibration  behavior  '(that  is,  when  the  modal  shape  possesses 
a  small  number  of  axial  waves),  the  resulting  motion  in  the  first  approxi¬ 
mation  for  "freely  supported"  ends  is  again  nonperiodic,  reflecting,  in 
the  present  development,  an  energy  level  equal  to  that  for  the  periodic 
solution.  Although  the  solution  corresponding  to  a  small  number  of 
axial  waves  was  not  the  general  case  under  investigation,  it  nay  be 
conjectured  that  the  introduction  of  a  more  general  deflection  function 
than  the  one  employed  will  result  in  a  nonperiodic  motion  reflecting  a 
lower  energy  level  than  the  periodic  solution.  For  the  representative 
cases  investigated  herein,  the  nonperiodic  motions  for  an  assumed  modal 
shape  possessing  either  a  large  or  small  number  of  axial  waves  appear  to 
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be  almost  periodic.  This  phenomenon  has  extreme  significance  in  the 
interpretation  of  experimental  results.  Unfortunately,  very  little 
experimental  information  concerning  the  nonlinear  free  vibration  be¬ 
havior  of  thin  shells  exists. 

When  the  number  of  circumferential  waves  is  small,  the  Karman-Donnell 
theory  becomes  invalid  and  it  is  necessary  to  employ  a  more  accurate 
set  of  strain-displacement  relations,  specifically  those  attributed  to 
Sanders,  while  at  the  same  time  retaining  the  midsurface  inertia  terms. 

Due  to  the  inherent  complexity  of  this  nonlinear  shell  vibration  prob¬ 
lem,  the  solution  presented  herein  is  based  only  upon  the  basic  chess¬ 
board  modal  pattern.  This  then  results  in  a  vibration  behavior  which 
exhibits  pronounced  hardening  characteristics.  Retention  of  additional 
terms  in  the  expression  for  the  radial  deflection  function  will  permit 
the  possibility  of  the  occurrence  of  a  nonperiodic  motion  which  may,  for 
certain  representative  cases,  appear  to  be  almost  periodic,  exhibiting 
either  a  hardening  or  a  softening  characteristic. 

For  the  case  when  a  large  number  of  circumferential  waves  occurs  (Karman- 
Donnell  theory),  and  for  the  case  when  a  small  number  of  circumferential 
waves  occurs  (Sanders  theory),  limitations  on  the  two  theories  are  pre¬ 
sented  based  upon  the  magnitude  of  the  allowable  rotations  at  large  dis¬ 
placements.  For  the  free  vibration  problem,  relatively  large  amplitudes 
can  be  tolerated  without  exceeding  limitations  on  rotation  magnitudes. 

Future  effort  on  the  problem  of  the  vibrating  circular  cylindrical  shell 
should  include  the  effects  of  damping  and  should  consider  the  problem  of 
the  nonlinear  response  of  a  shell  to  forced  vibration.  However,  the  re¬ 
sults  of  a  theoretical  development  are  most  beneficial  when  they  can  be 
complementary  to  the  results  of  am  experimental  investigation.  The  future 
development,  as  regards  the  behavior  of  the  nonlinear  free  and  forced 
vibration  of  a  thin,  circular  cylindrical  shell,  will  rely  heavily  upon 
the  findings  obtained  from  experimental  studies.  Certainly,  at  the  least, 
the  predicted  phenomenon  of  nonperiodic  free  v'bration  behavior  for  the 
circular  cylindrical  shells  should  be  verified  in  the  laboratojy. 
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Figure  1.  Circular  Cylindrical  Shell  Coordinate  System  and  Sign  Convention, 
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Figure  3.  Variation  of  Average  Maximum  Displacement  With  Frequency  Ratio 

for  the  Fixed  Parameter  Solution  (Reference  16,  u  =  0.50  and  2.00). 
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Figure  U.  Variation  of  Average  Maximum  Displacement  With  Frequency  Ratio 
for  the  Fixed  Parameter  Solution  (Reference  16,  rj  =  1.0). 
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Figure  5.  Comparison  of  the  Maximum  Radial  Displacements  for  the  Fixed  Parameter  Solution  (Reference  i6) 
and  the  Present  Free  Parameter  Solution  as  a  Function  of  the  Time  Parameter  t 


Figure  6.  Maximum  Radial  Displacement  for  the  Present  Free  Parameter  Solution  as  a  Function  of  the 
Time  Parameter  t  for  the  Minimum  Energy  Case. 


Figure  7.  Variation  cf  Strain  Energy  Parameter  With  y,  for  the 


Present  Free  Parameter  Solution  (y,  --  1.0) 
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Figure  9.  Variation  of  Strain  Energy  Parameter  With  y.,  for  the 
Present  Free  Parameter  Solution  (y^  =  0.01)  . 
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"  igure  10.  Effect  of  Finite  Deflections  on  Frequency  Ratio  Using  the  Sa^de^s 
Strain-Displacement  Relations  (p  =  1.0) 


Figure  11.  Effect  of  Finite  Deflections  or.  Frequer.r.-  Ra- 
Strain-Displacement  Rela*.  l.r.s  ~  C.  v 
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Figure  lU.  Modal  Lines  for  a  Full  Wave  Aspect  Ratio  of  Unity. 
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Figure  15.  Variation  of  Maximum  Rotation  Magnitudes  With  Thickness-to-Radius 
Ratio  for  Varying  Values  of  Wave  Amplitudes  (all  u  and  t)  =  1.0)  . 


Figure  16.  Variation  of  Maximum  Rotation  Magnitudes  With  Thickness-to-Radius 

Ratio  for  Varying  Values  of  Wave  Amplitudes  (p  =  0.50  and  tj  =  1.00) 
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Figure  17.  Variation  of  Maximum  Rotation  Magnitudes  With  Thickness-to-Radius 
Ratio  for  Varying  Values  of  Wave  Amplitudes  (p  =1.0  and  t)  =  1.0) 


Figure  18.  Variation  of  Maximum  Rotation  Magnitudes  With  Thickness-to-Radius 
Ratio  for  Varying  Values  of  Wave  Amplitudes  (p  =2.0  and  tj  =  1.0) 
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Figure  20.  Variation  of  Maximum  Rotation  Magnitudes  With  Thickness-to-Radius 
Ratio  for  Varying  Values  of  Wave  Amplitudes  (p  =  1.00  and  t]  =  0.0V 
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APPENDIX  I 


DEVELOPMENT  OF  THE  APPROXIMATE  SOLUTION  FOR  A  LARGE  NUMBER  OF' 
CIRCUMFERENTIAL  WAVES  FOR  THE  FIXED  PARAMETER  CASE 


This  appendix  presents  a  step-by-step  development  leading;  to  the  solu¬ 
tion  presented  and  discussed  in  earlier  sections. 
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The  governing  differential  equations  were  developed  earlier  and  are 
listed  again  as 
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The  deflection  function  for  a  shell  vibrating  in  the  nth  mode  is  given 

by 
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Substitution  of  the  deflected  shape  into  Equation  (58)  and  subsequent 
solution  yield  the  following  expression  for  the  stress  function: 
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The  deflection  function  and  the  stress  function  given  by  Equations  (59) 
and  (60)  lead  to  satisfaction  of  the  requirement  that  the  circumferential 
displacement  be  a  periodic  function  of  the  circumferential  coordinate. 
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The  Lagrangian  (L  *  T-U)  can  be  written  as  follows: 
l  2nR 
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The  relations  between  the  stresses  and  the  stress  function  are  given  by 
Equation  (48).  Substitution  of  the  deflected  shape  and  the  stress  func¬ 
tion  into  Equation  (6l)  for  the  Lagrangian,  after  integration  over  x 
and  y  ,  yields 
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where  t  is  a  nondimensional  independent  variable  related  to  the  time 
t  by  the  following  relation: 
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The  application  of  the  variational  principle  relative  to  y^(x)  yields 
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the  following  Euler  equation: 
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Equation  (64)  is  a  single  second-order  differential  equation  in  the  in¬ 
dependent  variable  x  .  By  making  a  change  of  variable,  the  single 
second-order  differential  equation  in  t  is  reduced  to  two  coupled 
first-order  differential,  equations  in  x  .  The  required  change  of 
variable  is  given  by 
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Then,  the  two  coupled  first-order  equations  in  x  become 
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Solution  of  the  two  coupled  differential  equations  was  carried  cut  by 
standard  numerical  techniques  in  conjunction  with  utiliration  o.’  "he 
Burroughs  5500  computer  in  the  Stanford  University  Computation  Center. 


51 


APPENDIX  II 


DEVELOPMENT  OF  THE  APPROXIMATE  SOLUTION  FOR  A  LARGE  NUMLER  OF 
CIRCUMFERENTIAL  .WAVES  FOR  THE  FREE  PARAMETER  CASE 


This  appendix  presents  a  step-by-step  development  of  the  solution  when 
the  assumed  deflected  shape  is  more  general  than  that  used  in  Appendix  I. 


The  deflection  function  is  given  by 
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The  stress  function,  obtained  as  in  Appendix  I,  is  given  by 
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The  continuity  condition  for  the  circumferential  displacement  is  given 
by  the  following  mathematical  representation: 
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Substitution  of  the  deflected  shape  and  the  stress  function  into  the 
continuity  condition  yields  the  following  relation: 
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Hence,  the  deflected  shape  now  becomes 
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Substitution  of  the  deflected  shape  and  the  stress  function  into  the 
Lagrangian  (in  a  manner  exactly  analogous  to  that  described  in  Appendix  1 
after  integration  over  x  and  y  yields  .  , 
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The  application  of  the  variational  principle,  relative  to  y^(x)  and 
y3(t)  ,  yields  two  Euler  equations.  After  the  introduction  of  two 

changes  of  variable  defined  ty 

^1 
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and 
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the  two  second-order  Euler  differential  equations  are  transformed  into 
four  first-order  differential  equations.  The  first  two  are  given  by 
Equations  (74)  and  (75).  The  remaining  two  equations  are 
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Solution  of  the  four  coupled  differential  equations  was  carried  out  by 
standard  numerical  techniques  in  conjunction  with  utilization  of  the 
Burroughs  5500  computer  in  the  Stanford  University  Computation  Center. 
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APPENDIX  III 


DEVELOPMENT  OF  THE  APPROXIMATE  SOLUTION  FOR  A  SMALL 
NUMBER  OF  CIRCUMFERENTIAL  WAVES 


This  appendix  presents  the  development  of  the  differential  equations  of 
motion  by  the  use  of  the  strain-displacement  relations  deduced  by  Sanders. 


The  deflection  functions  are  given  by 
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Substitution  of  the  deflected  shapes  into  the  strain  energy  and  the 
kinetic  energy  expressions  (Equations  (17)  and  (16),  respectively)  after 
integration  over  x  and  y  yields  the  following  expression  for  the 

Lagrangian: 
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The  application  of  the  variational  principle  relative  to  y  (t)  ,  y ._  ( t )  , 
>  y-j('r)  ,  y9(t)  ,  y^T)  ,  and  y^(T)  yields  seven  Euler  equation 
The  following  changes  of  variables  are  now  employed: 
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Next,  the  seven  second-order  differential  equations  (in  the  independent 
variable  t  )  are  transformed  to  fourteen  first-order  differential 
equations.  The  first  seven  of  the  fourteen  equations  are  given  by 
Equation  (80).  The  remaining  seven  equations  are 
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and 


Solution  of  the  lU  coupled  differential  equations  was  carried  out  by 
standard  numerical  techniques  in  conjunction  with  utilization  of  the 
Burroughs  5500  computer  in  the  Stanford  University  Computation  Center 
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